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The present paper is concerned with diffusion processes running 
on tubular domains with conditions on nonreaching the boundary, re- 
\ spectively, reflecting at the boundary, and corresponding processes in 

the limit where the thin tubular domains are shrinking to graphs. The 
methods we use are probabilistic ones. For shrinking, we use big poten- 
tials, respectively, reflection on the boundary of tubes. We show that 
there exists a unique limit process, and we characterize the limit pro- 
cess by a second-order differential generator acting on functions defined 
(-H , on the limit graph, with Kirchhoff boundary conditions at the vertices. 

> 

1. Introduction. The present paper is concerned with diffusion processes 
running on tubuiar domains with Dirichlet (i.e., absorbing-like) (resp., Neu- 
mann, i.e., reflecting) boundary conditions, and the respective processes ob- 
tained in the limit where the thin tubular domains shrink to graphs. Prob- 
I lems of this type have been intensively studied before in the case of Neumann 

■rj" ■ boundary conditions, both by probabilistic tools [21, 22] and analytic tools 

'nI" I [2, 8-10, 12, 13, 15, 38, 41]. The case of Dirichlet boundary conditions was 

known to present special difficulties, which explains why there have been, 
! up to now, fewer works concerned with this case, and, in fact, these are only 

■ concerned with either special graphs or special shrinking procedures, lead- 

ing mainly (with the exception of [2, 9, 10, 12]) to limiting processes which 
"decouple at vertices" [7, 11, 15]. 

Before explaining these difficulties and entering into details let us motivate 
the reasons to undertake such studies, pointing out also some connections 



■ with other problems and giving some historical remarks. 
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In many problems of analysis and probability one encounters differen- 
tial operators defined on structures which have small dimensions in one or 
more directions. Let us mention as examples the modeling of fluid motion 
in narrow tubes, or in nearly two-dimensional domains (see, e.g., [42]), the 
propagation of electric signals along nearly one-dimensional neurons (see, 
e.g., [3, 7, 11]), the propagation of electromagnetic waves in wave guides [31], 
the propagation of quantum mechanical effects in thin wires (in the context 
of nanotechnology) ; see, for example, [2, 9, 10, 12, 13, 15, 17, 24, 32, 33, 35, 
41, 48]. Such geometrical structures tend in a certain limit (mathematically 
well described in general through a Gromov topology) to a graph. Modeling 
dynamical systems or processes on such structures by corresponding ones 
on a graph might present certain advantages (e.g., PDEs becoming ODEs 
on graphs; more dimensional spectral problems reduced to one-dimensional 
ones). In any case the study of dynamics and processes on graphs can be 
considered as an idealization or a "first approximation" for the study of the 
corresponding objects in more realistic situations. 

There is a rich literature on differential operators on graphs. Diffusion 
operators and evolution equations were considered originally in work by 
Lumer [37], and subsequently by many authors; see, for example, [5, 40, 
49, 50]. Elliptic and parabolic nonlinear equations on graphs have been dis- 
cussed, for example, in relations to applications in biology, for example, 
in [11]; see also, for example, [3, 7] for nonlinear diffusions on graphs in 
connection with neurobiology. Heat kernels on graphs have been studied 
in particular in [39]. Hyperbolic nonlinear equations on graphs have been 
studied, for example, in [31]. 

In quantum mechanics, Schrodinger equations on graphs are considered 
as models of nanostructures; see, for example, [6, 17, 32, 33]. Work has 
been particularly intense in the study of spectral properties of Schodinger- 
type operators on graphs; see, for example, [24, 32, 33, 35]. Such models of 
quantum mechanics on graphs also play an important role in the study of the 
relation between classical chaos and quantum chaos; see, for example, [16, 
24, 35, 43, 44]. 

For the study of the limit of differential operators on thin domains of 
(and corresponding PDEs) degenerating into geometric graphs (and corre- 
sponding ODEs) we refer to [30, 42, 50] and especially to the surveys by 
Raugel [42] (which discuss topics like spectral properties, asymptotics and 
attractors). For the study of parabolic equations and associated semi-groups 
and diffusion processes we also refer to [42]. Corresponding hyperbolic prob- 
lems in connection with the modeling of ferroelectric materials have been 
discussed, for example, in [1]. 

Probabilistic methods for the study of processes on thin domains of M" 
have been developed by Freidlin and Wentzell in the case of Neumann bound- 
ary conditions. They exploit the consideration of slow, respectively fast, com- 
ponents going back to [20], applied to the thin tubes problem [21]. In these 
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studies the basic probabilistic observation is that for a Brownian motion in 
a thin tube along a line, the component in the transverse direction is fast, 
and the one in the longitudinal direction is slow. The control in the limit 
exploits the assumption on the reflecting properties of the fast component, 
together with a projection technique onto the longitudinal direction. In [21] 
it is shown that the diffusion coefficient for this limit process is obtained by 
averaging the diffusion coefficient for the process in tubular domains with 
respect to the invariant measure of the fast component with suitable changed 
space and time scales. 

Analytically the Laplacian in the transverse direction has a constant 
eigenvalue (ground state in the transverse direction), which then yields 
a natural identification of the subspace of L? — over the thin tube cor- 
responding to the eigenvalue for the Laplacian in the transverse direc- 
tion with the — space along an edge. Results about this approximation 
concern convergence of eigenvalues, eigenfunctions, resolvents and semi- 
groups [13, 15, 25, 38]. Besides, operatorial and variational methods also 
methods of Dirichlet form theory have been used [8] . 

The identification stressed above is no longer possible in the case of Dirich- 
let boundary conditions on the boundary of the thin tube, since the lowest 
eigenvalue of the Laplacian in the transverse direction diverges like 1/e^, 
where e > is the width of the narrow tube. (For a probabilistic study of 
the first-order asymptotics of the lowest eigenvalue of the Dirichlet Lapla- 
cian in tubular neighborhoods of submanifolds of Riemannian manifolds, 
see [28].) This has been pointed out clearly and posed as an open problem 
by Exner (see [4]). In order to nevertheless manage analytically the limit to 
a graph, one has to perform a renormalization procedure, first introduced 
in [2], and extended in [9, 10], for the case of a V-graph (waveguide). More 
general cases with Dirichlet boundary conditions have been managed in the 
case where the shrinking at vertices is quicker than the one at the edges; 
however, then one has "no communication between the different edges" (i.e., 
"decoupling") on the graphs; see [25, 38, 41]. The interest in discussing the 
case of Dirichlet-boundary conditions is particularly clear in the physics of 
conductors, where such boundary conditions arise most naturally, both in 
classical and quantum mechanical problems. However, in the other type of 
applications we have mentioned there is also an interest in studying bound- 
ary conditions that are different from the Neumann ones, since boundary 
conditions influence the limit behavior, and one is interested to obtain on the 
graphs the most general possible boundary conditions at the vertices (even in 
the case of an "A^-spider graph" there are AT^-different possible self-adjoint 
realizations of a Laplacian on the spider; see, for example, [17, 29]). 

The present paper mainly discusses the case of shrinking by potentials, 
and the goal is to determine the limit process on a given graph. This shrink- 
ing by potentials corresponds to confining the process in thin tubes around 
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the graph, not reaching the boundary almost surely, and in this sense is re- 
lated with Dirichlet boundary conditions (the latter property corresponding 
however to a completely absorbing boundary). In Sections 2 and 3 we con- 
sider special cases, because the consideration of these cases illustrate better 
the methods we use. 

In Section 2 the case of a thin tube Jl^ in M" shrinking to a curve 7 in M" 
is discussed. The tube has a uniform width e > 0. In the tube we have 
a nondegenerate diffusion process X"^ with a drift consisting of two parts, 
one continuous and bounded, the other of gradient type, pushing away from 
the boundary, so that the first hitting time of at the boundary dVf^ is 
infinite almost surely. We also construct a diffusion process X on 7 and 
show (Theorem 2.2) that if ^'"(0) converges weakly to ^(0), then also X^ 
converges weakly to X. If pathwise uniqueness holds both for X^ and X, 
then X^ also converges to X almost surely as e 4- 0- We also state correspond- 
ing results for a process in VL^ with a reflecting boundary condition on the 
boundary dVf^ (Theorem 2.3). These results are obtained in a similar way as 
those obtained by our shrinking with potentials in the first part of Section 2. 

In Section 3 we discuss the case of shrinking thin tubes in to an 
A'^-spider graph in M". In this section, we often use the methods discovered 
by Freidlin and Wentzell [21], extend their method to the case of diffusion 
processes instead of Brownian motions and apply it to the case of shrinking 
by potentials. The process X^ in the domain Vf" consisting of N tubes is 
defined in a similar way as in Section 2, e > being the parameter of shrink- 
ing to the A^-spider graph F for e 1 0. We prove again that the first hitting 
time of X'^ at the boundary dO.^ is infinite and that the laws of {X^ : e > 0} 
are tight in the topology of probability measures on C([0,+oo)), if their 
initial distributions are tight. We then show that any limit process is strong 
Markov and study the transition probabilities from the vertex O to any edge 
of the spider graph F. This requires quite detailed estimates of the behav- 
ior of the process X^ in a neighborhood of O in Q,^ . These results imply 
that the boundary condition at O should be a weighted Kirchhoff bound- 
ary condition for the functions in the domain of the generator of the limit 
processes X. (This is one of the types of boundary conditions known from 
the general discussions on boundary conditions for processes on graphs; see, 
for example, [12, 17, 29, 32-34].) The weights are determined explicitly from 
the construction, as transition probabilities to the edges (Lemma 3.7). This 
is crucial to determine the generator of the unique limit process X (Theo- 
rem 3.8). Similar considerations lead to corresponding results for the case 
where X^ is a diffusion in Q!^ with reflecting boundary conditions on dVf 
(Theorem 3.9). 

In Section 4 we state the results in the case of thin tubes around gen- 
eral graphs, which are obtained immediately from the results in Sections 2 
and 3. These are systems consisting of thin tubes around finitely ramified 
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graphs in M" with edges which consist of C^-curves. Theorem 4.1 presents 
a result similar to the one for an A^-spider graph, showing, in particular, 
convergence of the diffusion process X'^ not leaving the system i}^ of tubes 
around the general graph to a diffusion process X on the graph. Again its 
generator is determined and an extension is given to the case of a diffusion 
with reflecting boundary conditions on dil.^ . Since the latter result is not 
only for a Brownian motion in the thin tubes, but also for reflecting diffu- 
sion processes in the thin tubes, it is also an extension of previous results of 
Freidhn and Wentzell [21]. 

All random variables discussed in the present paper are defined on a prob- 
ability space with probability measure P, and E[-] denotes their expectation 
with respect to P. For a locally compact topological subspace A of M^, let 
CoiA) := {/ G C{A) : lim|,|^+^ f{x) = 0}. 

2. The case of curves. In this section, we consider shrinking of thin tubes 
to curves. Let n be an integer larger than or equal to 2. Let 7 G C'^(M;]R") 
such that I7I = 1 [with 7 the derivatives of t — )• 7(t), and | • | the norm 
in M"], and assume that 7 has no self-crossing point, and 7 is a bounded 
function with a compact support. Let e > 0, (•, •) be the inner product on M", 
and d{x,'y) be the distance between x and 7. Note that d{x,^) is Lipschitz 
continuous in x. Define domains {i^^} by 

Jl^:={2;G]R":d(x,7)<e}. 

Consider a differentiable function u on [0,1) such that 

14(0)= 0, u'>0 limn'fi?) = +00 and -lim:; — = +00. 

Rn Rn iog(i - R) 

For example, if we define n(r) := r"/(l — r") for r G [0, 1) where a > 0, then u 
satisfies the conditions above. Let 

U'{x) = u{e-U{x,-f)), xen'. 

For e > 0, consider a diffusion process given by the following equation: 



/■t/\(,- rt/\(, 
X'it) = X'{0)+ a{X'{s))dW{s)+ h{X'{s)) 

(2.1) 



ds 



{VU'){X^{s))ds, 



where -'^^(O) is an f2^-valued random variable, W is an n-dimensional Wiener 
process, a G Cb(M''; M"), b G Cfe(M";M") and C is the first hitting time 
of at the boundary dQ^ of Q^. Let a := aa"^ (with a'^ the transpose 
of o"), and assume that a is a uniformly positive definite matrix. Then, the 
solution of (2.1) exists uniquely; see, for example, [47]. 



Lemma 2.1. = +00 almost surely for small e > 0. 
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Proof. Assume n > 3. Note that does not hit 7 almost surely in 
this case. Let be the solution of (2.1) replacing X'^ifi) and by x 
and C|> respectively, where Cx is the first hitting time of X^ at It is 
sufficient to show that C| = +00 almost surely for x near to dVt^ . By the 
tubular neighborhood theorem and Theorem 1 in [18], there exists a C^- 
diffeomorphism (/> = {4>i,4>2) from f]^ \ 7 to {y = (^1,^2) G K x :0 < 

\y2\ < e} which satisfies, for small e, 

= 7""*" o 7r(x) and (f)2{x) = d{x , j)V d{x , , xGf]^\7, 

where 7r(x) is the nearest point in 7 from x. Note that (/> is a C^-function on 
Q'' and (Vvr, V?/^) = for small e. Hence, {S/<l)i,VU^) = and V02Vf/^ = 
e~^n'(e~^(i(-,7))V(i(-,7). By Ito's formula, we have 



(2.2) + / VMXUsmXUs))ds 

Jo 

+ 2E / a,,(X|(s))5,a,0i(X^(s))ds, 
MXm) = M^)+ / V02(^.^(s))cT(X|(s))dVF(^.) 

JO 

(2.3) ^° 

+ 2E / 'aii(^l(s))5.9,</'2(X^(5))ds 

-e-i / u'{e-'d{X',{s),j))Vd{;j)\x.is)ds. 
Jo 

Moreover, again by Ito's formula, 

\Mxm)\' 



02 



»|2 + 2 / {UXl{s)),VUXl{s))a{Xl{s))dW{s)) 



+ 2 / {MX'As)),VMXUsmXUs)))ds 




ftACJ / " \ 



DIFFUSION PROCESSES IN THIN TUBES 7 

-2e-i / \(p2{Xl{s))\u'{e~^d{Xl{s),^))ds 
Jo 

+ / tr^ce[VUXl{s))a{Xl[smMXl[s))a{Xl{s))f] ds. 
Jo 

Let 

d := sup{\{VMxMx)f^\^:xe G {y € : |y| = 1}}, 

b:= sup ^2((/)2(2;), V(/)2(x)6(x)) + l^2{x), ^ aij{x) didj(j)2{x)^ 

+ trace[V(^2(x)a(x)(V(^2(x)f7(x))^]^ . 

Take cq G (0, 1) such that sup^gj^^ — 2xu'{x)) < and 

f(x):= \ exp -2 / ^-^dz\dy, x^%e^). 

Then, by Ito's formula, for 5 such that < 5 < 1 — cq and for x such that 
co£ < (^(a^jT) < ^(1 ~ ^)-, we have that 

i?[/(|02(X^(T^°^AT-(i-^)))|2)]</(d(x,7)'), 
where T"" := inf{t > : 7) = c} for c> 0. Since 

£;[/(|</>2(x^(r^°^Ar-(i"^)))|2)] 

= /(cge2)P(T^o^ < T^(i-^)) + /(e2(i _ 5)2)p(r^o£ > ^£(1-5)) 

and 

we have 



/(e2(l-5)2)_/(c2,2)- 

The assumptions on u imply that f{£^{l — (5)^) diverges to +00 as 5—7-0. 
Hence, the proof is achieved from the fact that converges to C| as 

In the case where n = 2, since X^ can hit 7, we need a little arrangement. 
Let and fil be the two domains consisting of fi^ \ 7, and 6'^{x) be 1 if 
X S ri^, —1 if X G and if a; G 7. By the tubular neighborhood theorem 
and Theorem 1 in [18] again, there exists a C^-diffeomorphism 0= {(l)i,4'2) 
from il^ to {y = (1/1,^2) G M x (— e,e)} which satisfies, for small e, 

(f)i{x) = o 7t(x) and (f)2ix) = 0^ {x)d{x,j), x G 0^, 
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such that (2.2) and (2.3) hold. Thus, we can discuss this case in a similar 
way as the case where n > 3. □ 

Theorem 2.2. Define a diffusion process X by the solution of the fol- 
lowing equation: 

X{t) = X{0) + A o ^~\X{s)){j o j-\X{s)),a{X{s)) dW{s)) 
Jo 

+ A o 7-i(X(5))(7 o 7-i(X(s)), 6(X(s))) ds 
Jo 

(2.4) + ^ ° r\Xis))\aiX{s)f^ o ^-\Xis))\'ds 

+ A o ^-\X{s)){a{X{s)f^ o j-\X{s)), 
Jo 

a{X{s)f^o^-\Xis)))ds. 

Note that X is uniquely determined as a process on 7. 

If X^{0) converges to a ^-valued random variable X[{)) weakly, then the 
process X^ converges weakly to X in the sense of their laws on C([0, +00); M*^) 
as e 1 0. 

Moreover, if pathwise uniqueness holds for (2.4) and (2.1) for all £>Q, 
and -'^^(O) converges to a ^-valued random variable X{{)) almost surely, 
then X^ converges to X almost surely, as eJ^O. 

Proof. Note that equation (2.2) holds even if we replace x and C| 
by X^, X^[{)) and C^, respectively. Lemma 2.1 implies 

(2.5) sup d(X^(t),7) ^0, eiO, 

<G[0,+oo) 

almost surely. Hence, the boundedness of the coefficients implies the tight- 
ness of the process (j)i{X'^). Let X be any limit process of subsequence of X^ . 
Then, we have X S C([0, +00); 7) almost surely by (2.5). Hence, taking e 1 
in (2.2) with replacing X|, x and C| by X^, X^(0) and C"^, respectively, 

^i{X{t)) = + f V0i(X(s))a(X(s)) dW{s) 

Jo 

+ / V^i{X{s))h{X{s))ds 
Jo 

+ 2E / aij{X{s))d,d,(t>i{X{s))ds, 
where W is an Wiener process. 
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Noting that (pi{X{-)) is a stochastic process on M and \ V(j)i{x)a{x)\ > for 
rr E 7, the law of (/>i(X(-)) is uniquely determined by this equation; see The- 
orem 3.3 of Chapter IV in [27]. Applying Ito's formula to 7((/)i(X(t))) and 
noting that 7((^i(X(-))) = X{-), dicpi = 7j o 7"^ on 7 for z = 1, 2, . . . , and 
didj(j)i = (7i 07-1 )(7j 07-1 ) + (7j 07-1 )(7j 07-1 ) on 7 for i,j = 1,2,...,N, 
we have that X satisfies (2.4); therefore, the first assertion holds. The second 
assertion is obtained in a similar way. □ 

The argument above is also available in the case where the boundary dfi'^ 
carries a Neumann boundary condition, for the generator of the process, in 
the following sense. Consider a diffusion process X^ which is associated with 

in and reflecting on dJl^. Then, X^ can be expressed by the following 
equation: 

(2.6) X^{t) = X'{{))+ f a{X^{s))dW{s)+ [\{X' {s)) ds + <^%X'){t), 

Jo Jo 

where is a singular drift which forces the reflecting boundary condition 
on (9i7^; see [46]. Discussing this case in a similar way as above, we obtain 
the following theorem. 

Theorem 2.3. Define a diffusion process X by the solution of the fol- 
lowing equation: 

X{t) = X{0) + A o 7-i(A(s))(7 o j~\X{s)),a{X{s)) dW{s)) 
Jo 

+ r 7 o 7~'(A(s))(7 o ^~\X{s)),b{X{s))) ds 
Jo 

(2.7) ^11^"^° j-\X{s))\aiXis)f^ o j-\X{s))\'ds 

+ r 7 o ^-\X{s)){a{Xis)f^ o r\Xis)), 
Jo 

a{X{s)f^oj-\x{s)))ds. 

If X^{0) converges to a j-valued random variable X{0) weakly, then the 
process A^ converges weakly to X in the sense of their laws on C([0, +oo); M*^) 
as e 4, 0. Moreover, if pathwise uniqueness holds for (2.7) and (2.6) for 
all e > 0, and A'^(O) converges to a '^-valued random variable X{Q) almost 
surely, then X^ converges to X almost surely, as e|0. 
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Remark 2.4. In this section, the shape of tubes was taken to be cyhn- 
drical and the "confining" potential has been defined by the scahng of 
a fixed function U . However, neither the shape of the tubes nor the scaUng 
property are essential. If is "along 7" (in the sense that the gradient of 
is normal to the tangent of 7), the same results hold. In the case where 
is not along 7, some effect of C/^ remains in the limit process; see [19, 45]. 

3. The case of AT-spiders. In this section, we consider the shrinking of 
thin tubes to A^-spider graphs. The argument in this section is the main part 
of this article. Consider an n-dimensional Euclidean space M", let d{-,-) be 
the distance function in R" and let O be the origin. Let {ci^f^-^ be N different 
unit vectors in M" and li := {scj : s € [0, 00)}. Consider an A^-spider graph T 
defined by F := Ui^i -^i- T is also called an A^-star graph. Let A be the set 
in M" given by 

A:= [j {xGM":(x,ei) = (2;,ej)}. 

For X G M" \ A, let tt{x) be the nearest point in T from x. Note that tt{x) is 
uniquely determined for all x € M" \ A. 

Let Ui be given similarly to u in Section 2 for i = 1,2, . . . ,N (so that Uj de- 
termines the potential acting in the thin tube around /,). Let Cj be a positive 
number for i = 1,2, . . . , N , 



K :=max{2\/2cj/Y'l - {ei,ej):i,j = 1,2, N,i ^ j} 

and kq G {0,k,). Ci has the interpretation of width of the tube around /j. 
Let U he a, function on M" with values in [0,oo], and assume 

U{x) = Ui{c^'^d{x,T)), X G {x G M"':7r(x) G Ii,d{x,Ii) < a, \x\ > k}, 

U{x) = +00, X G {x G M" :7r(x) G Ii,d{x,Ii) > q, |.x| > k}, 

U{x)<+oo, xG{xGM":|x| <Ko}, 

Vt := {x:U{x) < 00} is a simply connected and unbounded domain, dO, is 
a C^-manifold and U\q is a -function in Q,. This structure Q is sometimes 
called a "fattened" A^-spider. In addition, we assume 

U(Xm) 

lira {-VU{xm),Vd{xm,dQ)) = +00 and - lim - — — - — ^^i— — = +00 

for any sequence {xm} which converges to a point x G d^l. Define domains 
{ni:i = l,2,...,N} in by 

ili := {x £ Q \ A : tt{x) G !{, \x\ > k} 

for i = 1, 2, . . . , A^. Let := eQ, Qf := eQi, and U'{x) = U{e-^x) for x G M" 
for ah e>0. Note that U^{x) G [0, +00) for x G 0^, dU^ is a C^-manifold, 
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and U^\Qe is a C^-function on Q,^ . Consider a diffusion process X'^ given by 
the following equation: 

X'{t) = X%0)+ / a{X'is))dWis)+ / b{X'{s))ds 

(3.1) 

(VC/^)(X^(s))ds, 

where -'^^(O) is an fi'^-valued random variable, is the first hitting time 
of at 9f]^, W is an n-dimensional Wiener process, a S Cb(M";M" ^ M") 
and 6 G Cfe(M";]R"). Define a stochastic process X^ by the solution of (3.1) 
with replacing X^{0) by x, and P| by the law of X^ on C([0, oo); M"). Let 
a(x) := (T(x)(T'^(rc), and assume that a is a uniformly positive definite matrix. 
Define a second-order elliptic differential operator L on il.^ by 

172 r\ r\ Tl r\ 

then the generator of X"" is a closed extension of (L — S/U^ ■ V) in L^(r2^, dx) 
for any e > 0. Since a is a uniformly positive definite matrix, the process X^ 
exists uniquely for all e > 0. 

The following lemma implies that X^ does not exit from Vt^ almost surely. 



Lemma 3.1. = +oo almost surely for all e > 0. 



Proof. Locally, the discussion in the proof of Lemma 2.1 is available. 
Hence, by using the strong Markov property of X"^, we have the assertion. 
□ 



Next we shall study the tightness of {X^ :e > 0}. 

Lemma 3.2. // the laws of {X^ (0) : e > 0} are tight, then the laws of 
{X^ :e > 0} are also tight in the sense of laws on C([0,oo);M'"). 

Proof. In view of Theorem 2.1 in [21] it is sufficient to show that for 
any p > there exists a positive constant Cp such that for all y € M" there 
exists a function fp on M" which satisfies the following: 

(i) fjliy) = I, fjiix) = Ofor\x-y\>p and < f^ < 1. 

(ii) {fp{X^{t)) + Cpt;t > 0) is a submartingale for sufficiently small e. 

Now we choose fp and Cp satisfying the conditions above. Fix p > 0, and 
take eo > such that Eq < p/{1Qk). When y G fi^o (where denotes the 
closure of in R") and |?/| > /9/2, choose fji G C°°(M") such that: 
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• /^(x) = /^(7r(x)) for X G \ A and /^(x) = for \x-y\> p/A; 

• /p^(?/) = 1, < /p^ < 1, ||V/||oo < 8/p and HWlloo < 64/p2. 

Since fp{x) = for |x| < 2Keo and V7r(2;)VJ7'^(x) = for |x| > 2Keo, it follows 
by Ito's formula that 

fyix%t))- f Lfy{x^{s))ds 

Jo 

is a martingale for all e < Eq. Hence, choosing Cp larger than (8//9 + 64//>^) x 
(||o"||^/2 + llfelloo), conditions (i) and (ii) are satisfied for e < eo- 
When y G and \y\ < p/2, choose E C°°{W) such that: 

• fy(x) = Pp{-K{x)) for X G \ A, fy{x) = 1 for |x| < p/4, and /^(x) = 
for |x — y| > p; 

• /p^(?/) = hO<fp< 1, ||V/||oo < 8/p and || WlU < 64/p2. 

Here, note that 2Ke < /5/4 for e < Eq. Similarly to the case where y £ 
and \y\ > p/2, one proves that conditions (i) and (ii) are satisfied for e < Eq 
with the same Cp as above. 

When y ^ IW, choose G C°°{W) such that /^(y) = 1, /^(x) = for x G 
0*^0, and /p satisfies condition (i) above. Since moves in $7^, fp{X^{t)) = 
for all t and e < Eq- 

Thus, for all p > 0, {fp : y £ M"} and Cp are chosen in such a way that 
conditions (i) and (ii) are satisfied. □ 

Now, we assume the tightness of {X^{0) :e > 0}. By Lemma 3.2 we can 
choose a subsequence {X^ : e' > 0} of {X^ : e > 0} such that the laws of its 
members converge weakly in the sense of laws on C([0, oo); M*^). Define X as 
the limit process of this subsequence, and to simplify the notation denote the 
subsequence e' by e again. Prom now on we fix X as the limit process of X^ . 

For tj; G C([0,+oo);M"), let f%w) := inf{t > : |t(;(t)| = c} and T^^^w) := 
mf{t>0:w{t)^A,\7T{w{t))\ = c} for c> 0. 

Theorem 2.2 determines the behavior of X on F \ O. Hence, to charac- 
terize X , we need to determine the boundary condition for X at O. Now 
we give some lemmas. The following lemma implies that the edge which X 
goes to, starting from O, is independent of the edge which X comes from. 
Therefore, we obtain in particular that X is a strong Markov process on F. 

Lemma 3.3. Let {5{e) : e > 0} be positive numbers satisfying the condi- 
tion that limeioE~'^5{E) = +oo. For B G ^{W) denoting the Borel 
subsets ofW^], 

sup{|P^^(u;(T^(^)) G S) - P5(w;(T^(^)) G S)| :x G |x| < 3Ke} 
converges to as e^O. 
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Proof. Define a process X| by the solution of the equation 

ft ft 



(3.2) 



Xl{t) = x+f a{eXl{s))dW{s) + e f h{eXl{s))ds 
Jo Jo 

- f\vU){Xl{s))ds 
Jo 



for x G $7 and e > 0, where W is an n-dimensional Wiener process defined 
by W{t) = £~^W{£'^t) for t £ [0, oo). It is easy to see that the law of {X^{t) : 

t > 0) is equal to {£-^XI^{e'^t) : t > 0) for x G Letting P| be the law of X| 
on C([0,oo);M"), we have 

(3.3) Pxi^it) G dx) = PI^{£~^w{E'^t) £ dx) 

for t G [0, oo), x G O and e > 0. By (3.3), it is sufficient to show that 

(3.4) |P^^('u;(r^(^)/e) e £-1^) _ p^[w{T^^'^l') G £-^B)\ ^ 

as £ tends to 0, uniformly in xG {yG ^'■\y\ < 3k}. Define stopping times 

To{w) := mi{t > : w{t) ^ A, |7r(u;(t))| > 3k}, 

fkiw) := mi{t > Tk-i:w{t) ^ A, |7r(u;(t))| > 4ft;}, /c G N, 

Tk{w) := M{t > fk:w{t) ^ A, \'K{w{t))\ < 3ft;}, k G N, 

for we C([0,oo);M"). Note that |u;(rfe)| = 3ft; for A; = 0, 1, 2, . . . , and \w{fk)\ = 
Ak for /c = 1,2,3, .. . almost surely under for x G and |a;| < 2ft;e. Since 
A7r(x) = and V7r(x)VC/(x) = for |x| > 2ft;, Ito's formula implies 

n{Xl{t)) = n{fu{Xl))+ f ^ Vn{Xl{s))a{£X^{s))dW{s) 

(3.5) 

+ e f ^ V-n{Xl{s))b{£X^{s))ds 

Jrk(X-) 

for t G [ffc(X|),rfc(X|)], x eVL and |x| < 3ft;e. Since the diffusion coefficient 
of the one-dimensional process |-7r(X|(t))| is uniformly elliptic, and T^^'^^^ 
diverges to infinity as e J, almost surely under , there exists a sequence 
{77(e)} converging to as e 4, such that 

sup P^"(r^(")/" <r3«) <7?(e). 

\x\=Ak 

On the other hand, since aa^ is uniformly positive definite, X| hits {x G : 
|x| < 5'} with positive probability for all x G e > 0, 5' > 0. Hence, let- 
ting a(e) be a sequence of positive numbers such that Q(e) < 2ft;, and a(e) 
converges to as e 4,0, we obtain that 

p{£):= inf p^^(r"(^) <r4'^)>o 

|x|=3k 
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for all e > 0, and that p{e) converges to as e J, 0. Moreover, we have 

oo 

fc=l 

oo ^ ^ 

=e/ / 

-'{xieC: |7r(xi)|=3K} J{j/i6n: |7r(j/i)|=4K} 



X • • • X P^^^('u;(r4'') G dyi,T^'' < r"(")) 



<r?(.)J](l-p(e))^ 
k=l 

p{e) 

Hence, if rj{e)/p{e) converges to as e 1 0, P^(r'5(")/^ < r"(")) converges to 
as e 1 0. Now we choose a^e) so that 77(e) /p{£) converges to as e | 0. Then 
pe{T&{e)/e ^ fa{e)^ converges to as e 10. Thus, for (3.4), it is sufficient to 
prove that 

(3.6) sup |P^^('u;(r^(^)/^)Ge-^S)-P5('u;(r^(^)/e)G £-^5)1^0 

\x\<a{e) 

as e 4, 0. To show this convergence, we use the coupling method. Let a' G 
Cb°°(M'";M" (g)M") and 6' G Cfe°°(R";M") for / = 1, 2, . . . , such that 

lim sup \a\x) — a{x)\=0 and lim sup \b\x) — b{x)\ = for Af > 0. 

Let X be fixed, and consider a pair of stochastic processes (X|'',Xq') defined 
by 

X^J{t) = x+ f a\eXlHs))dW{s) 
Jo 

+ e [ b\eX^/{s))ds- [ {VU){X^/{s)) ds, 
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+ e / h\eX''J{s))ds- 

2a\ex2)-^{xi - X2){xi - X2f {a\ex2)-^Y 



iVU)iX's\^))ds, 



where 



H'^\xi,X2):=In 



\a'-{ex2) ^{xi-X2)\^ 



for xi,X2 € M", and /„ is the unit matrix. Note that {Xx\Xq'') is uniquely 
determined because of the smoothness of and . We define 

V{y) := \y\~'^y, ^>^''(xi,X2) := a^{£Xi) - a^{ex2)H^'^{xi,X2), 
<f^\xi,X2) := eh\exi) - {VU){xi) - eb\ex2) + (V[/)(x2) 

for y G € • ^ 2ko} and xi,X2 €z {z G M" : \z\ < kq}- Similarly to the 
argument in Section 3 in [36], there exists a positive constant K such 
that {xi, X2)^V {xi - X2)P} > K for xi,X2 € {z G M" : \z\ < kq} for 

small e. By the equi-continuity of {a'}, we can choose p G (0, 2ko) satisfying 

2(xi - X2, ^'^''(xi , X2)) + trace($^''(a;i, X2)$^''(xi, X2)^) 

- \<^>'^\xi,X2fV{xi - X2)P < K/12 

for jxi — X2I < /O, |xi| < Kq, \x2\ < Ko,l = 1,2,3, . . . (see [36]). For e' € [0,p), 
define a stopping time on C([0,oo);R") x C([0,oo);M") by 

^''{w,w') :=mi{t>0:\w{t) -w'{t)\ ^ {e,p), 

1^(1)1 > Kq, or > Kq} 

for w,w' £ C([0,oo);M"). By Ito's formula and the choice of p, we have 

p'/'P{\X^J{:^^\x^/,X^S'))-^o\^''i^^^^^ 

<i^[|X,^''(=^^'(X^'',X^''))-X^''(^^'(X^'',X^''))|2/3] 



\Xl;\s) - X^o\^) 
^\¥'\Xl:\s),rs\^)f 



-4/3 



3 



3--' {Xl\x'ii'') 



\Xl:\s)-X^S 



el, 



UO 



-4/3 
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y<V{Xl\s)-X^A\smds 



el, 



- ' ' 18 



\X^J{s)-X^S\^)\~'''ds 



^ El^^'{X^J,X^J)]- 



18//3 



Hence, letting e' 4, 0, we have the following two estimates: 
(3.7) 

18/5^/3 



(3.8) E[^\X^J,X's')]< 
On the other hand, by Ito's formula. 



K 



E 



\x^/{^"{x^/,x^S'))-^\-Y 



E 







\Xl\s)-x\-^] \Xl\s) - x\-%,^. 



2 / , '""^ -{\X%'\s)-x\>Ko/2} 

X <J 2{Xl\s) - x,eh\eXl>\s)) - VU{Xl\s))) 



+ trace [ex' (eX^'' (s))cj' (eX|'' (s)) 



I , N ^a;' (s) X 



<y\eXl'\s)) 



\Xl'\s)-x\ 



{|XS''(s)-x|>«;o/2} 



(s) - j;| 



<CE[^\Xl^\X'J)], 



where z+ := max{0, z} for 2; G R and C is a positive constant independent 
of / and X. This inequality together with (3.8) implies 

72Cp4/3 



(3.9) 



P{\Xl\^^\Xl\X^J))-x\ = K,)< 



kIK 



\x 



|2/3_ 
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Similarly, we have 

(3.10) P{\rs\^\Xl\rs'))\ = Ko) < ^-^^\x\'/', 

where C is a positive constant. Noting that and Xq" converge to X| 

and Xq in law as I — > +00, respectively, for each e, by the coupling inequality 
(see [36]) we have 

sup \P^{w{T^^''^'') G e^^B) - PS{w{T^^'^/') G e^^B)\ 

\x\<a{e) 

< sup sup\P{X'J{T^^''>/')ee~^B)-P{X'J{T^^'y')ee-^B)\ 

|x|<a(e) I 

< sup supP(|X^''(^°(X,^'',X^''))-X^'^(=^0(X^'M^''))|/0) 

|a;|<a(e) I 

< sup supP(|X^''(^°(Xi'',X^''))-X^'^(^0(X^'M^''))|=p) 

kl<a(£) ' 

+ sup supP{\X^J{^'^{X'J,X'o^))-x\ = Ko) 

\x\<a{e) I 

+ sup supP(|X^''(=^0(l^'M^''))| = Ko). 

\x\<a{e) I 

This inequality, together with (3.7), (3.9) and (3.10) yields (3.6). □ 
The next lemma implies that O is not absorbing for X. 
Lemma 3.4. 

f E[I{,.,\,\^S'}{X{s))]ds = 0{6') 
Jo 

as 6' iO, for aUt>0. 

Proof. To simphfy the notation, let X^{0) = e 0^. It is sufficient to 
show that 

E[I^,.,\^^,)\^S'}{X{s))]ds = 0{6') 

as 5' 1 0. By Fatou's lemma, we have 

ft 

E[I{x:\jT{x)\<5'}{X{s))]ds 

(3.11) <liminf /*E[I{,,|,(,)|<3,,}(X^(s))]ds 

+ liminf / E[^^,3^,<\^(^)\<s,y{X%s))]ds. 
£4.0 Jo 
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To show that the second term is 0{6') as 5' J, 0, let / be a continuous 
function on M such that I{xm:3Ke<x<5'} <f< ^xm:2Ke<x<25'} and F{x) : = 
lo fo f^^) dzdy. Noting that tt{x) = (cj, x)ei for x S and i = 1, 2, . . . , A'", 
we have V7r(x)7r(a;) = 7r(x) for x ^Vt such that |x| > 2k. Since V7r(a;)VC/''(x) = 
and A7r(x) = for x ^Vt^ such that \x\ > 2k£, we have 



E[Fi\7riX%t))\)]-F{\7rix^)\) 



1 
2 

+ 



E 



E 



f{\<X^is))\) 



\<X-{s))\ 



ds 



ds. 



It is easy to see that is dominated uniformly in e > 0. Moreover, 

it holds that < F' < 26' and < F{x) < 26' x for x G M+. Thus, by uniform 
ellipticity of a = aa^ , we have the following estimate: 

y{\7TiX^ism]ds<C6' 



for some constant C. Hence, 

ft 

E[I{x:3Ke<Hx)\<5'}iX'{s))] ds = 0{6') 



lim inf 

£4-0 Jq 



as 6' 1 0. This yields that the second term of (3.11) is equal to 0{6') as 6' J, 0. 
The proof is finished by showing that 

(3.12) / i?[I{,,|,(,)|<3,,|(X^(s))]ds = 0(e) 

Jo 

as e|0. Define stopping times {T|,f|} by 

:= mf{u > T^_i{w) : \tt{w{u))\ > Ake}, k€N, 

:= inf{n > f^{w) : \tt{w{u))\ < Sks}, A; G N, 

for 'u;gC([0,oo);M"). Then, 
ft 

E[l{x:\^{x)\<3Ke}iX'is))]ds 

f^Kl T'''%w)P^{dw))p^4w{rD e dx,4<t) 



< 



< sup 

x£{yen:\7T{y)\=3Ke} 



[T'^%w)P^{dw))f2PxM<t) 
^ k=l 
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By using the notation in the proof of Lemma 3.3, we have 

sup I T^^'{w)P^{dw) = e^ sup f T^^{w)P^{dw). 

x£{yGQ : \n{y)\=3Ke} J xG{j/Gf2 : |7r(y)|=3K} J 

It is easy to see that 

sup sup / T'^''{w)P^{dw) < +00. 

e>0 xG{ygn : |7r{j/)|=3K} J 

Hence, for (3.12), it is sufficient to show that 

(3.13) 5]P,^.(f|<t)<Ce"i 

k=l 

for some constant C. For w S C([0, 00); 0), let Mt{w) be the number of transi- 
tions of w from the set {x : |'7r(x)| = 3k} to the set {x : |'7r(x)| = 4k} 
during the time interval [0, t] . Then, 

00 „ 

(3.14) Y.P-^^^k<t) = / K~H{w)P^-rAdw). 
k=i 

Take / G C°°([0,oo)) such that / > 0, < /' < 1, f" > 0, supp/" C 
[2k, 3k], f{x) = for X < 2k and /(3k) < /(4k). Define Y^'' by 

for 2; G and i = 1, 2, . . . , A^. Since {ei,VU {x)) = for x G {VL-i : \x\ > 2Ke}, 
by Ito's formula we have 

Y^'\t) = f{{e,,x)InAx)) 

+ [ /'((e„X^(s))In.(X^(s)))(e„c7(eX^(s))cil?(5)) 

+ e [' f\{e,J^,{s))In,{X^AsmeiMeXUs)))ds 



1 



+ f"{{e,,Xl{s))lnAXl{s)))\a{eXl{s)fe,\''ds. 
It is clear that 

TV 

E[M,-2,{Xl^,^^)\ < J2 sup E[K-.,{Y^'% 

j^—l x:\n{x)\<4:H 

where Mtiw) is the number of up-crossing of w for the interval [/(3k), /(4k)] 
during the time interval [0, i]. Hence, by (3.13) and (3.14), it is sufficient to 
show that 

(3.15) sup E[K-2t{Y,'n]<Ce~' 

x:\tt{x)\<Ak 
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with a constant C for alH = 1, 2, . . . , A^. Let i be fixed and m G N. Define 
and ffc by 

TO := 0, 

TO ■.= mf{u> 0:Y^''{u) < f{3K)}, 
fk := mf{u > rfc„i : Z^'^n) > 
Tk:=mf{u>fk:y^'\u)<f{3K)}, 

Then, 



A; EN. 



k=l 
m—1 

k=l 

m 

Y,E[Y^\fkAt)-Y^'\Tk-iAt)] 
k=l 

m—1 



k=i 

m—1 



k=l 

Since /" > 0, we have 



TkAt 
TkAt 



J2E[n'\fk/\t)-n'\n-iAt)] 



k=l 



(3.16) 

< E[Y^'\t A fm)] + E[Y^''\t A To)] + Ciet 
with a positive constant Ci. Let 

f* := maxjffc : < A; = 1, 2, 3, . . .}, 
n := max{Tfc : < t, A; = 1, 2, 3, . . .}, 

M{t) := r/'((e.,X,^(s))%(X,^(s)))(e„a(eX^(5))dW^(s)). 



When n < Y^'\t) < f{4K). When n < n, /"((e„ X|(s))]Io, (X|(s))) 
for s G [f^,t]. Thus, we have 

Y^'\t) = Y^'\n) + Mit)-M{n) 



Hence, 
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+ e f f'{{e,,Xl{s))lu, {Xl{s))){e,MeK{s))) ds. 

J r* 



Y^^\t) < /(4k) + 2 sup |M(s)| + C2et 

0<s<t 



for \x\ < 4/t with a constant C2. By the Bur kholder-D avis-Gundy inequahty 
we have 

e\ sup y/'Xs)l < /(4k) + 2C3\/t + C2et 

for |x| < 4k with a constant C3. Thus, letting m— )• +00 on (3.16), we have 
for |x| < 4k 

(/(4k) - f{3K))E[AftiY^'')\ < 2/(4k) + AC^Vt + (Ci + 2C2)et. 
Therefore, replacing t by e~^t, (3.15) is obtained. □ 

The lemmas above yield that the boundary condition at O is a weighted 
Kirchhoff boundary condition. Hence, the next step is to determine the 
weights associated with the edges. Let Y^ be a diffusion process defined 
by the solution of the following stochastic differential equation: 

(3.17) Y^{t) = x + a{0)W{t)- f {VW){Y^{s))ds. 

Jo 

Note that is a special case of with the condition X^{0) = x, and Y^ 
does not hit almost surely. Denote the law of Y^ on C([0,oo);M") by Qf^. 
It is easy to see that the law of Y^ is the same as that of £Y^_i^ (e~^) • By (3.2) 

one has that the law of X| converges to that of Y^ as e | 0, and therefore, 
the law of X^ and that of Y^ are getting closer as e | 0. In particular, we 
have 

]im\PS{w{T^) G ni) - Qh{w{T^) e nt)\ = 

£4,0 

for all and i = l,2,...,A^. Since this holds for all c > 0, it is possible to 
choose a subsequence of e (denote the subsequence by e again) and positive 
numbers /3(e) which satisfy lim£^o/3(e) = +00, and 

(3.18) lirn|P5(w;(r^(-)^) G ni) - Qh{w{T^^'^n G f^DI = 
£4.0 

for i = 1,2, . . . , N . Let 6{e) := e/3(e). Then, 6{e) satisfies the conditions in 
Lemma 3.3. 

Now we assume that cr(0) = In where In means the unit matrix. This 
assumption enables us to determine the weights of the edges explicitly. Let 



Pi 
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We remark that when Uj is independent of i, then we have pi := ^ / 
(l2iLi hence the weights {pi} are determined by the ratio of the area 

of the cross-section around the edge /j. Then, the foUowing lemma holds. 

Lemma 3.5. If a{0) = In, then 

lim sup |P^^(u;(r'^(")) e!^f) -pil =0 

fori = l,...,N. 

Proof. Applying Lemma 3.3 to both and X^, and using (3.18), it 
is sufficient to show that 



(3.19) 



\im\Qh{w{T'('^)eni)-p,\=0 



for i = l,...,N. 

We make a similar discussion as in the proof of Theorem 6.1 in [21]. 
Let be the invariant measure of the Markov chain {Y^^rf.)}, where r| are 
stopping times defined by 



0, 



inf{n > Tf,_^{w) : \7t{w{u))\ > 6{e)}, keN, 

inf{n > f^{w) : \tt{w{u))\ < Sks}, /c G N. 

Define a measure /.t^ on by 

fi'^{dx) ■.= exp{—W{x))dx, xGil.'^, 

a function space by {/ G C^(i7^) : lim^,. ^(j, ^qe-j^q f{x) = 0} and a bi- 

linear form by 



{Vf{x),Vg{x))p'{dx), /,5G^('^^ 



Then, the pre-Dirichlet form {S"^ , ^^^S"^)) on L'^{Q^,p^) is closable, and 
is associated to the Dirichlet form obtained by closing {S"^ , Note 
that is an invariant measure of Y^; see [23]. By Theorem 2.1 in [26] we 
have 

^f{B)= / v^{dx) / / lB{w{t))dt Qlidw) 

for B G ^{W). Let Bf := {x G nf:6{e) < \tt{x)\ < 26{e)}. Then, 



(3.20) 



{x£Q° : |7r(x)|=3«;e} 



u^{dx) 



{w{t))dt 



QUdw) 
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u'{dx) I In,{w{T^^'^))QUdw) 

On the other hand, let 

Z{t) := -6{e) + W{t), f := inf{t > : \Z{t)\ > 25(e) - Sks}, 
where Vl^ is a one-dimensional Wiener process starting from 0, and 

/X j-y 
-25(e) J-2S(e) 

Then, by Ito's formula we have 



{dw) 



E[F{Z{f))]-F{-5{e)) = \E 



'5{e),5(e)]{Zt)dt 



Since F can be computed explicitly, we see that F{—5{£)) = and 
E[F{Z{f))] = F{25{e) - ?>ne)P{Z{f) = 25{e) - S^e) 
6{e) — ?,n£ 



A5{e) - 6Ke 



[25{eY + 25{e){5{£) -'ine)]. 



Thus, it follows that 

e' 







~S(e),S(£)]{Z{t))dt 



■25{ef + o{5{ef). 



On the other hand, the strong Markov property and the reflection principle 
imply that 



{w{t))dt]Ql{dw) = E 



~S{e),S{e)]{Z{'t))dt 



for all y G {x € fi^ : |'/r(a::)| = (5(e)}, because the left-hand side is independent 
of the behavior of w moving in G il'" : |7r(a;)| > 5(e)} under Qy. Hence, it 
holds that 

(3.21) / (/^ lBi{w{t))dt^Ql{dw) = 25{ef + o{5{ef) 

for all y G {x G : |vr(x)| = 5(e)}. By Lemma 3.3, (3.20) and (3.21), we have 
^x^{Bl) = {25{ef + o{5{ef))u^{{x G fi^' : |^(x)| = 3Ke}) 

(3.22) 

x(Q-oKr'^'^)ef^f) + Oe(i)). 

Since Q}){w{T^^^^) G 17f ) = 1, we have, as e i 

1 ^ 

(3.23) v^{{x G : |7r(x)| = ^ne}) = -5(e)-2 ^/.^(Sf ) + o,(l). 

1=1 
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Dividing both sides of (3.22) by those of (3.23), we obtain that 
QUw{T'^'^) G in) = f + o,(l). 

By the definition of fi^, the continuity of a and b, and cr{0) = In, ^"^(^f 
can be expressed explicitly as 



where cOn-2 is the area of the (n — 2)-dimensional unit sphere. Therefore, (3.19) 
is proved. □ 



The statement in Lemma 3.5 can be improved as follows. 
Lemma 3.6. 



hmlim sup \P^iw{T^')ent)-p,\=0 
S'lO eiO |a;|<3K£ 



for i = 1, . . . , N . 



Proof. In view of Lemma 3.3, it is sufficient to show 

limlim|P5(u;(r^')GOf)-p,|=0 
o'iO eiO 

for i = 1, 2, . . . , A^. Define stopping times {t|, f|} by 
T^{w) :=0, 

f^(w) := inf{n > tI_i{w) : \-k{w{u))\ > 6{e)}, keN, 
rliw) := inf{n > f^{w) : \tt{w{u))\ < Ske}, A; G N. 
By the strong Markov property, we have 

oo „ 

(3.24) =Y, ^^TS'y{w)PS{dw) 

k=l-' 

X j P^-(r^'<r3-)fe(y)P^(.._^)(u;(r'^(^))Gdy) 

and 



k=l 

oo „ 

(3.25) =^J I^^^_^^^s'y{w)PS{dw) 



fc=i' 



X I p,Pl{T'' < r3-)P^(,._j(u;(r'5(^)) e dy) 
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for i = 1, 2, . . . , A^. Let /il and be functions on [0, oo) given by 

, , , , 2{ei,b(x)) 
'i_(^):=max sup ^ — ^, 

* a:6Qf:|7r{a;)|=max{^,3K£} Fl^J 6i| 



:= min 



inf 



2(ei,6(a;)) 



- xeQf :|7r(x')|=max{^,3Ke} |cr(x)'^eiP ' 

respectively. Define functions si and on [0,oo) by 



ht{z")dz" dz' 



h%{z")dz"j dz', 
respectively. Then, for y £ {x : |vr(x)| = 5{e)} we have 



st{\7r{w{T^' A r3^^))|)P„^((iu;) - st{6{e)) 



sli{e,,w{T'' AT^^')))P^idw) - si((e„y)) 






/ii(w'(s))k(i(;(s))^e, 

f>«)(s) 



T„ |2 
« I 



+ 



X exp 



(ei,6(u;(s))) exp 



/o 



ht{z')dz' ] ds 

w{s) 



ht{z)dz'\ ds 



pe 

y 



<0. 

Hence, it holds that 



for y £ {x gQ'^ : |'/r(x)| = 6{e)}. Since 



we have 
(3.26) 



for y £ {x : |7r(x)| = (5(e)}. Similarly we have 



(3.27) 



s%{5{e)) - s%{3Ke) 
s%{5') - s%{3Ke) 
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for y G {x G $7^ : |vr(2;)| = 5{e)}. Let Mrps' {Xq) be the number of transitions 
of Xq from the set {x G il^ : |vr(x)| = Ske} to the set {x €fi^ : \tt{x)\ = S{e)} 
during the time interval [0,T^' (X^))]. By Lemma 3.5, (3.24), (3.25), (3.26) 
and (3.27), we have 

si((^(e)) -si(3Ke) 



< 



si(5')-si(3Ke) 

oo „ 



oo „ 

I — 1 ^ 



k=l 

s%i6{e)) - s%{3Ke 
~ s%i6') - s%i3Ke) 

+ °^(^) .1(5')-.! (3k.) ^^^tAXo)]- 
By the definitions of si and s^, we obtain 

,.,(aUm)-AA^ s%(S{e)) - s%{3k8) \ 

'™r*<^' [ .i(.--)-»i(3..) - J-""- 

where C is a constant independent of 6' , and for each 5' > 

si (<^(£)) -.1(3^5) 

.i(^0-^i(3^e) 

On the other hand, a similar discussion as in the proof of Lemma 3.4 implies 

E[M^s'iXh)] = Sie)-'os^il). 

Therefore, we have 

limsuplimsup(P5('u;(r'^') G nf) - pi) < 0. 

S'iO eiO 

Similarly we obtain 

limsuplimsup(pi - PS{w{T^') G Of)) < 0. 

5'iO eiO 

These inequalities yield the conclusion. □ 



We need a little more improvement of Lemma 3.6 as follows. 
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Lemma 3.7. 



limlimlim sup I P^fwfT^ ) g Of) =0 

S'iO SiO eiO 1^.1 <5 



fori = l,...,N. 

Proof. In view of Lemma 3.6 it is sufficient to show 

limlimlim sup \P^{w{T^') eni)-pi\=0. 

5'iO SIO elO 3f,s<\x\<5 

By Lemma 3.6 again, 
\P^{w{T'')eQI)-p^\ 

[ Py{w{T^') G l^f)P^"(w;(r3«=) G dy,T^^' < T^') 

i{ye!l=: |7r(i/)|=3Ke} 

<Pi\P',{T'^''' <T^')-1\+ P^iT^""' > T^')Ini{x) + o,,s'{l). 

Here, o^^s'i^) means a term which converges to as 5' ^ after letting e ^ 0. 
Hence, it is sufficient to show for 6' > and i = 1,2, ... ,N 

(3.28) limlim inf P^JT^""^ <T^') = 1. 

510 elO xeni:3Ke<\x\<S 

Let T^{w) := mf{t > 0:w{t) = O} and fix i. By Theorem 2.2 the law of 
(T3«^(X|.),r^'(^lO) converges to that of {T0{X^),T^' (X^)) as e i for 
x"^ G {y G il^ -.Ske < \y\ < 6} such that converges to x G li, where the 
process Xx is determined by the following stochastic differential equation: 

Xx{t)=x+ [\e„a{Xx{s))dW{s))+ [\e„b{Xx{s))) ds, 
Jo Jo 

te[0,TO{Xx)AT^'{Xx)]. 
By using Jj = ne'>o U£<£' ^^'^ compactness of {y G M" : \y\ <6}, we have 
lim inf PUT^'''<T^') 

eiO xGOf :3ree<|x|<<5 

= inf P{tO{X.,)<t''{Xx)). 

x£li:0<\x\<S 



Since crcr^ is uniformly positive definite, we have 
lim 

SiO x£ 

This proves (3.28). □ 



lim inf P{TO{X^)<T^'{Xx)) = 1. 

&iO x£U:0<\x\<6 
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The lemmas above determine the boundary condition for X at O. Now 
let us characterize X by a generator of a process on V. Let 

dej{x) := lim -{f{x + sei) - f{x)) 

s^O S 

for any differentiable function / on Ij and i = 1,2, . . . , N . Define a second- 
order differential operator £j on /j by 

(3.29) A := i|a^(x)e,|2a2 + {b{x),ei)de^ 

for i = 1,2, . . . , N . Define the second-order differential operator C on Co{T) 

:= | / G Co{r):f\j^\o e Ciil, \ O) for aU i = 1, 2, . . . , iV, 

\\m.Cif{sei) has a common value for i = 1, 2, . . . , A^, 



5^K(hm(aeJ)(se,))=0 , 

i=l ) 

Cf{x):=Cif{x), xeh\0, 
Cf{0) ■.= \nnCif{sei). 

Note that Cf{0) does not depend on the selection of i = 1, 2, . . . , A^. We 
call {pi} the weights of the Kirchhoff boundary condition at O, and call 
Z]i^iPj(l™s^o(c^ei/)('Sei)) =0 the weighted Kirchhoff boundary condition 
at O. 

Theorem 3.8. Consider diffusion processes defined by (3.1). As- 
sume that (j{0) = In and the law of X^{0) converges to a probability measure 
/io on T. Then, X^ converges weakly on C([0, +oo);M"') to the diffusion pro- 
cess X as e J, 0, where X is determined by the conditions that the law ofX(0) 
is equal to fiQ and 



(3.30) E 



f{X{t))-f{Xis))- / CfiX{u))du 







for t > s >0 and f € S'{C), where {^t) is the filtration generated by X. 
Therefore, C is the generator of X . 

Proof. From Lemma 3.2 we have that {X^} is tight. We are going to 
show that there is a unique limit point in this family. Let X be any limit 
point of {-'^^l, and denote the sequence converging to X by {X^} again. 
Since this martingale problem is well-posed (see [21]; [14] for the relationship 
between martingale problems and partial differential equations, and [37] for 
the uniqueness of the semigroup generated by £), it is sufficient to prove 
that X satisfies (3.30). Fix s > 0. Let 6' be a positive number. Define the 
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following stopping times: 



To 
To 
fk 
Tk 



mf{u>s:X{u)=0}, 
mf{u>Tk^i:\X{u)\>5'}, A; G N, 
mf{u>fk-X{u) = 0}, /cGN. 



Then, for fe&{C), s<t 



E 



E 



f{X{t))-f{X{s))- / £f{X{u))du 

J s 

oo ^ 

^(/(X(tAffc))-/(X(tArfc_i)) 



.fc=i 
oo 



Cf{X{u))du 



fc=0 



f{x{t^Tk))- f{x{t^fk)) 



tATk 



Cf{X{u))du 



tAfk 



Because of Theorem 2.2 the second sum vanishes. We estimate the first sum 
as follows: 

^(/(X(tAffc))-/(X(tATfc_i))- / li^{X{u))du 



E 



< 



E 



J2 {f{X{fk))-f{X{Tk-l))) 
t 



+ l|-C/||oo^ 



^x:\x\<S'}iX{u))du 



+ sup |/(x)-/(0)|. 

|x|<5' 



Clearly, the third term on the right-hand side converges to as 5' ^ 0. By 
Lemma 3.4 the second term on the right-hand side converges to as ^ J, 0. 
The first sum on the right-hand side is equal to 

oo N 

Y.Y.(f(^'^i'^-f(0))P{X{fk)eh,fk<t\^s 



(3.31) 



k=l 1=1 



oo N 



k=l i=l 

Let 6 G (0, 6') and let, for any e > 0: 

T^'^ := mi{u > s : \tt{X'{u))\ < 5}, 



ff := inf{u > rlti ■ k(X^(n))| > 5'}, A; G N, 
:= mi{u > ff : \tt{X'{u))\ < 5}, G N. 



30 



S. ALBEVERIO AND S. KUSUOKA 



The distributions of the pairs {X^ ,f^'^ ,t^'^) converge weakly to those of 
(X, ffc,Tfc) as 5|0 after e|0. Hence, by Lemma 3.7 we have 

P{X{fk)£li,fk<t\^s) 

= limlimP(X^(f,^'^)GOf,ff <t|^.) 

= nmlim I P^{w{T'')€ni)PiX^i4'!^)€dy,ff<tl^s) 

= {pi + OS'{l))P{fk<t\^s)- 

Note that X^fcli-f ('^fe < t\^s) is equal to the expectation of the number of 
transitions of X from the point O to the set {x G F : = 6'} during the time 
interval [s,t] [with respect to a general initial condition X(0)]. Approximat- 
ing that by the expectation of the number of transitions of X^ from the set 
{x : \tt{x)\ = 6} to the set {x £ : |vr(x)| = 5'} during the time interval 
[s,t], similarly as in the proof of Lemma 3.4 we obtain the estimate 



oo 



Ct 



k=l 



'Pifk<t\^s)< ^ 

with a positive constant Ct depending only on t. Hence, by (3.31) we have 

{f{XitAfk))-fiX{tATk-i))) 



E 



Ct 
< — 



N 



V5'lim/'(sei)Pi + o(<5') 



Since / € ^(>C), the right-hand side converges to as (5' J, 0. □ 



Similarly as in Section 2, the argument above is also available in the 
case where the boundary of fi^ carries a Neumann boundary condition. 
Consider a diffusion process X^ which is associated to L in and satisfies 
the reflecting boundary condition on d^l^ . Then, X^ can be expressed by 
the following equation: 

(3.32) X'{t)=X'{0)+ f a{X'{s))dW{s)+ [\{X' (s)) ds + ^'{X'){t), 

Jo Jo 

where is a singular drift which forces the process to be reflecting on dfi'^; 
see [46]. Note that X^ depends on fi^ but is independent of . Discussing 
this case in a similar way as we did in the case of Dirichlet boundary condi- 
tion we obtain the following theorem. Let 
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:= | / G Co(r) : f\j^\o e \ O) for alH = 1, 2, . . . , iV, 
lim£j/(sej) has a common value for i = 1, 2, . . . , A^, 

N 

^ Pi(\im{deJ){sei) 



i=l 
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Cfix):=Cifix), xeh\0, 
Cf{0):=lhjiQf{sei), 

where £j is given by (3.29). Note that Cf{0) does not depend on the selec- 
tion of i = 1,2, . . . , N . 

Theorem 3.9. Consider the diffusion processes defined by (3.32). 
Assume that cr(0) = /„ and the law of X'^{0) converges to a probability mea- 
sure fiQ on r. Then, {X^} converge weakly on C([0, +oo); M") to the diffu- 
sion process X as e -10, where X is determined by the conditions that the 
law of X{0) is equal to fiQ and 



E 



f(X{t))-f{X{s))- / CfiXiu))du 







for t > s >0 and f G ^(^), where {^t) is the filtration generated by X. 
Therefore, C is the generator of X . 

Remark 3.10. The weights {pi} of the case of Neumann boundary con- 
dition can be obtained from the wights {pi} discussed in Theorem 3.8 in the 
heuristic limit where the potential Uj around each edge takes only the value 
on [0,1) and +oo on [l,-|-oo). 

Remark 3.11. As mentioned in Remark 2.4, we can discuss similarly 
the case where the shapes of the tubes are not cylindrical. However, 

if is not defined by a scaling of a fixed function U , the weights of the 
weighted Kirchhoff boundary condition cannot be determined uniquely. To 
handle this more general case, we have to assume that satisfies some 
uniform bound. 

4. The case of general graphs. In this section we present results obtained 
by combining the results of Sections 2 and 3, and, in this way, we cover more 
general graphs. Let A be a finite or countable set, H be a subset of A x A, 
{V\ : A G A} be vertices in M", {Ex x' : (A, A') G H} be C^-curves with ends 
{Vx,Vx'} and G:=U(A,A')GH^A,A'. Denote A ~ A' if (A,A')gH. 



Kx := max 
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Let us denote the length of i?A,A' by |£'a,A'|- Define (jx^yis) : s G [0, |£'a,A'|]) 
as the arc- length parameterization of -Ea,A' with 7a,a'(0) = ^A- Assume that 
the number of {Vx : A G A} n {x E M" : |x| < M} is finite for all M > 0, |£;a,A' I 
is finite for all (A, A') € H and 

lim(7A,Ai(s),7A,A2(s)> < 1 

for all A ~ Ai and A ~ A2 such that Ai 7^ A2. Let ca,a' be a positive number 
for (A, A') G H, and let 

2^/2cA,Ai/Jl-lim(7A ,Ai(s),7a,A2('5)) : Ai, A2 G A 

such that A ~ Ai, A ~ A2, Ai 7^ A2I 

for A G A. Let 7r(x) be a point in G which is nearest to x G M". Assume 
that there exists a small eo > and positive numbers {ka} such that tt(x) is 
uniquely determined for all x G Ua~a'{^ ^ • ^(^' -^a,a') < ca.a'^j d{x, Vx) > 
Kx£ and d{x,Vxi) > nys} and for all e G (0,eo]) and that 7a,a'(s) = 
sufficiently small s for each (A, A') G H. 

Let ua,a' be given similarly to u in Section 2 for (A, A') G H. For e G (0, eo]; 
let f/"^ be a function on IR" with values in [0, +00], and assume 

U^ix) =ux,x'icx),,e~^d{x,Ex,x')), 

X G {x G M":7r(x) G Ex,x' ,d{x, Ex^y) < cx,ye,d{x,Vx) > Kx£,d{x,Vy) > Kye}, 

U%x) = +00, 

X G {x G M":7r(x) G Ex,y ,d{x, Ex,y) > cx,ye,d{x,Vx) > Kxe,d{x,Vy) > Kye}, 

il.^ := {x:U^{x) < 00} is a simply connected domain, dil'^ is an (n — 1)- 
dimensional C^-manifold embedded in and C/^|n= is a C^-function on Q,^ . 
In addition, we assume 

lim {-\/U{xm),yd{xm,dn^)) =+00 and - lim - — — I^i— — = +00 

m^oQ [og[d{Xm,0\l'^)) 

for any sequence {xm} which converges to a point x G di}^. 

Consider a diffusion process given by the following equation: 

X^{t) = X%0)+ [ a{X^s))dW{s)+ [ b{X^s))ds 
(4.1) ^ ^° ^° 

''\vU'){X'{s))ds, 



where A'^(O) is an fi'^-valued random variable, W is an n-dimensional Wiener 
process, a G Cf,(M"; M" M") and b G Cb(M'*;]R"). Let a := aa'^ , and assume 
that a is uniformly positive definite. Define a second-order elliptic differential 
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operator L on 0,^ by 

i,j=i ■' 1=1 

Then is associated with [L — (VC/^, V)). Similarly to Section 3, it holds 
that does not exit from Vl^ almost surely. Assume that (riVx) = (y\In for 
all A G A where ax > 0. 

For (A, A') G H, define a second-order differential operator £a,A' on Ex^xi by 

1 _^ (f 

■= 2l^(^)^^A.^' °^a,a'(^)P^(/°7a,a')(7a,a'(^)) 

+ [{bix),'yx,x' °lx^'i^)) 
+ (o-(2;)^7A,A' o7a,v(2;), 

c^(2;)^7a,a' o Ta^a' (^))] ^ (-^ ° 7a,a' ) (7a,a' (^) ) ' 

for X G i^A.A' and / G C^(£'a,A') where s is the parameter for the arc- length 
parametrization 7a,a'- Let 

cl:x'lir--'^M-nx,xir))dr 

Ea; A~A lo '^""^ ^M-Ux:xi^)) dr 

By using these notations, define the second-order differential operator C 
on Co(G) by 

&{C) := |/ G Co(G) : /b, e ^'(^A,A' \ {l^A, l^A'}) for A ~ A', 

for A G A,lim£A,A'/(7A,A'('5)) li^s a common value for A' : A ~ A', 

^ ^'A,A'lim(^(/o^,,„(^)))=OforAGA}, 

A':A'~A ^ ^ ^ 

:= £a,A'/(x), X G i?A,A', (A, A') G H, 
£/(yA):= lim £A,A'/(:r), AG A, 

where the limit x ^Vx is along Ex^x' ■ Note that CfiVx) does not depend 
on the selection of A'. 

Since by locality the behavior of diffusion processes associated with dif- 
ferential operators is determined in a given point by the behavior in its 
neighborhoods, we have the following theorem by Theorem 2.2 and 3.8. 

Theorem 4.1. Consider the diffusion process X^ defined by (4-1)- As- 
sume that the law of ^^(0) converges to a probability measure fiQ on G. 
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Then, {X''} converge weakly on C([0, +00); M") to the diffusion process X 
as e^O, where X determined by the conditions that the law of X{0) is equal 
to Ho and 

rt 



E 



f{X{t))-f{X{s))- [ Cf{X{u))du 

J s 







for t> s>0 and all f G ^{C-), where {^t) is the filtration generated by X . 
The operator C as defined above is thus the generator of X. 

Similarly as in Sections 2 and 3, our discussion is also available for the 
case where the boundary carries a Neumann boundary condition for the 
process. Consider a diffusion process X^ which is associated with L in Vt^ 
and reflecting on dVf^ [defined similarly as the process described by (3.32)]. 

Let 



-A, A' 



E 



n-l' 



&{C) := |/ G Co{G):f\E^^,\{y^y^,^ G CliE^^x, \ {^a.^A'}) for A~ A', 

for A G A,lim£A,A'/(7A,A'('S)) has a common value for A' : A ~ A', 



^ ^^A,vHrn(A(/o^,,„(^)))=OforAGA}, 

A':A'~A ^ ^ 



Cf{x) 
CfiVx) 



-A, A 



fix) 



x^Vx 



XG^A,A',(A,A') GH, 

AG A, 



where the limit a; — >• Va is along Ex x' ■ Then, we obtain the following theorem. 



Theorem 4.2. Consider the diffusion process X^ defined above. As- 
sume that the law of X'^{0) converges to hq. Then, {X'^} converge weakly on 
C([0, +oo);M"') to the diffusion process X as e 4,0, where X is determined 
by the conditions that the law of X{0) is equal to ^0 ond 

ft 



E 



f{X{t))-f{Xis))- / Cf{X{u))du 







for t > s > and f G ^(>C) where {^t) is the filtration generated by X . The 
operator C as defined above is thus the generator of X. 

Remark 4.3. As mentioned in Remarks 2.4 and 3.11, similar discussions 
can be given for the case where the shapes of the tubes are not cylindrical. In 
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the case where a = In, 6 = 0, and -Ea,A' are straight, the result of Theorem 4.2 
coincides with Theorem 6.1 in [21]. 
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